Serial No. 09/714,292 



REMARKS 

STATUS OF THE CLAIMS 

Claims 1-16 were pending and stand rejected. 

By this Amendment, claims 1-16 are amended and new claims 17-19 are added. 
Therefore, claims 1-19 are now presented for consideration. 

No new matter is presented in the foregoing specification, claim and drawing 
amendments and new claims, and accordingly, approval and entry of same are submitted to be 
proper and are respectfully solicited. 

DRAWING OBJECTIONS 

In the Office Action, at page 2, item 1 , the drawings are objected to for various reasons 
provided therein. 

It is submitted that the objection to the drawings is now overcome because: a) one of the 
reference characters "502" has been amended to -501 --; b) the specification has been 
amended to conform with FIG. 29 with respect to reference character "il"; and c) the specification 
has been amended to conform with FIG. 7 with respect to reference character "a1 ." 

Moreover, Applicants traverse the objection with respect to "b2" of FIG. 9 not being 
mentioned in the description, since "b2" of FIG. 9 is mentioned in the specification at page 35, 
line 12. 

Reconsideration is respectfully requested. 

SPECIFICATION OBJECTIONS 

In the Office Action, at page 2, item 2, the disclosure is objected to because of 
informalities therein and, furthermore, the Examiner requested a copy of the convolution 
algorithm described at page 30, lines 16-18 and page 36, lines 19-22 and page 38, lines 23-24 
of the specification. 

Reconsideration is respectfully requested. 
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The specification is amended to overcome the objections thereto. 

Moreover, a copy of the convolution algorithm, as requested by the Examiner, is 
enclosed, as Attachment 2. 

CLAIM REJECTIONS UNDER 35 U.S.C. §112, SECOND PARAGRAPH 

In the Office Action, at pages 3-4, item 4, claims 1-16 are rejected under 35 U.S.C. §112, 
second paragraph, as being indefinite. 

Claims 1-16 are amended to overcome the rejection thereto. 

Reconsideration is respectfully requested. 

CLAIM REJECTION UNDER 35 U.S.C. §102 (e) 

Claims 1-3 and 5-16 are rejected under 35 U.S.C. 102(e) as being anticipated by Ogura 
et al. (U.S. Patent No. 6,502,984 B2) hereafter referred to as Ogura 984. 

Claim 1 is directed to an image processing apparatus, and recites "an image processing 
condition storing section to store an image processing condition when the radiation image is 
subjected to the image processing in accordance with a type of the photography device and 
a part of the target when the radiation image is obtained ... an image processing section ... to 
subject the radiation image ... to the image processing in accordance with the image 
processing condition." The image processing condition is stored when the radiation image is 
subjected to the image processing. That is, the radiation image or radiation image data exists 
and is subject to the image processing, the image processing being in accordance with the 
image processing condition. Furthermore, since the image processing apparatus of the present 
invention is provided with an image processing condition storing section, an radiation image is 
optimally processed based on the stored image processing condition, irrespective of the type of 
photographic device that took the radiation image and the part of the target as the subject of 
radiation image. 

Ogura 984 Reference 

Ogura 984 discloses "to provide an X-ray photographing apparatus which determines the 
X-ray irradiation time by measuring the body thickness of an object before X-rav irradiation , and 
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controls the moving grid on the basis of the determined X-ray irradiation time information to 
remove the influences of shadow images of lead foils formed in an X-ray image." (See Ogura 
984 at column 5, lines 53-59.) Thus, in the Ogura 984 apparatus, image processing is 
performed prior to X-ray irradiation such that an X-ray image is not subjected to the image 
processing in accordance with the image processing condition (see claim 1). This means 
that the Ogura 984 image processing is different from "the image processing," as recited in claim 
1. 

Furthermore, Ogura 984 is silent with regard to "an image processing condition storing 
section to store an image processing condition" associated with a type of the photography 
device and a part of the target," because in the Ogura 984 apparatus the image processing 
condition is measured (e.g., by measuring the body thickness). Thus, an image processing 
condition storing section is not disclosed or suggested by Ogura 984. 

Accordingly, claim 1 patentably distinguish over the cited art and is submitted to be 
allowable. 

Claims 2, 15 and 16 

Claims 2, 15 and 16, for reason similar to those of claim 1 are submitted to also be 
allowable. 

Claims 3-14 

Claims 3-14, which depend from claim 2, are also submitted to be allowable for at least 
the same reasons as claim 2, as well as for the additional recitations therein. 

CLAIM REJECTION UNDER 35 U.S.C. §103(a) 

Claim 4 is rejected under 35 U.S.C. 103(a) as being unpatentable over the combination 
of Ogura 984 and Ogura (U.S. Patent No. 6,314,198 B1) hereafter referred to as Ogura 198. 

It is submitted that the additional reference of Ogura 198 does not overcome the 
deficiencies of Ogura 984 because Ogura 198 does not disclose or suggest, in the portions cited 
by the Examiner, anything related to the recitation in claim 1 of "an image processing condition 
storing section to store an image processing condition when the radiation image is subjected 
to the image processing in accordance with a type of the photography device and a part of 
the target when the radiation image is obtained." 
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NEW CLAIMS 17-19 

New claims 17-19 are provide to afford a varying scope of protection. 

New claim 17 recites "an image processing condition table, having indices of parts of a 
body and types of photography devices used to image medical images, to store and output an 
image processing condition; and an image processor to enhance the medical image based on 
the image processing condition stored in the image processing condition table," and is submitted 
to be allowable. 

New claims 18 and 19 recite storing an image processing condition in a table having 
indices of parts of a body and types of a photography devices used to image medical images; 
and subjecting the medical image to processing based on the stored image processing 
condition," and are submitted to be allowable. 

Entry and consideration of this claim is respectfully requested. 

CONCLUSION 

There being no further outstanding objections or rejections, it is submitted that the 
application is in condition for allowance. An early action to that effect is courteously solicited. 

Finally, if there are any formal matters remaining after this response, the Examiner is 
requested to telephone the undersigned to attend to these matters. 

If there are any additional fees associated with filing of this Amendment, please charge 
the same to our Deposit Account No. 19-3935. Respectfully submitted, 



STAAS & HALSEY LLP 



Date: 



By: 



1201 New York Avenue, NW, Suite 700 
Washington, D.C. 20005 
Telephone: (202)434-1500 
Facsimile: (202)434-1501 



: 5^^^ 
Eric Berkowitz 
Registration No. 44,030 




I hereby certify that this correspondence is being deposit- 
ed with the United States Postal Service as first class 

— 
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ATTACHMENT 1 

TO THE AMENDMENT 
FILED APRIL 20, 2004 
IN SERIAL NO. 09/714,292 
DOCKET NO. 1080.1084 



ATTACHMENT 2 

TO THE AMENDMENT 
FILED APRIL 20, 2004 
IN SERIAL NO. 09/714,292 
DOCKET NO. 1080.1084 



©ISO/IEC 



FILTERING OPERATOR 



, • . ISO/IEC 12087,-2:1994(E) 
. ,convo!ye_2d 

COWOLVE.2D 

Description; . '>.•" . !' s <;.*.,:/:-..•» ; • ;v: y :., 

Performs twc-dimeniiobil convolution with a general impulse iesjttirise 'ftxacricm ajniy for several image 
boundary condition*. . ..... 



Name: 
convo!ve_2d 

Parameters; 

SRC 
DST 
SOB 

convolution 



destmafioh^jagite"" . v.- . / V • . 
impulse response array object ' 

convolution option \ 
upper left comer justified 

* J^^lxeli^'oexu^br"" 1 ' 
key pixel, reflected exterior 



ID 

ID 
ID 
SP 



.■.«.'. vl' 

:> \.; 
is-;: 



Objects: 

SRC structure: any image 

DST structure: , . m saxnc as SRC structure 

D$T dari type; '* 4 "saxnb'as SRC'dati tff£ 
SOB data type: SD.RD " ; ,l * 

SOB dimension: 2D 



Element functiooaHty: 

Index assignment: yes 

Match point: yes 

ROl conordfc yes 

ROI processing; no 

Source promotion: ' no 

Heterogeneous bands: no 

Chain able element: yes 

Asynchronous: yes 



. ISO)TEC.12087-2:19?4(E) © lsonEC 

. convoive_2d 

•Remarks:' ... 

1, The F1KS data object repository contains a set of fixed size impulse response arrays. 

2. Variable size impulse response arrays era be generated by £IKS took. 

. 3, .GneHdimensional convolution can be achieved with this operator by specifying an impulse response 
array with only a single row or column 

4. The impulse response array data type shall be SD if the source image dam type is ND or SD. 

5. The impulse response array data type shall be KD if the source image data rype is KD. 

Nomenclature; • 

SRC(x f y, z, t, b) Xs*Y$*ZgxT s xB s source image 

DST(x,y.z,tb) Xj^Y^xZ^Tj^B d destination image 

S(j. k, 1. m. n) J^xKjxLjxMsscNjj operator input image , tn * 

D(j*kJ,in.a) J^xKj^xL^xMdxNd operator pntpu i image ' " 

HQ;k)* J H xK H impulse response, funcdoa array . ... : . • . 

Qfe»^c^ loty pixel coordinate v •.- . 

The key pixel cwdinate is referenced to die engm'of '^ 'impuiiie response function array. It is act 
restricted to be within the bounds of the impulse respond favticn array. Le.. j k can be negative or 
larger than J H . ' ' - * t 

Definition: 

Upper left corner justified mode: 

The 2D impulse response function OTay..is..ro^^l^Oldegrees and translated with respect to' the >k 
pixel plane of the sour^OT£g£ A destination p«el is compute^or^li-spatid intersections of thb array 
and the aourc^imSge, The desdnarion image j-k pixel plane is* indexed froths upper left corner: The 
destmadpnr*image j-k plane is larger than the source image j-k plane:' ". N \ ; * 

Far/all ], m, zl and for \ . . 



The summation limits -axe 
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... COnvolvvJZd 

The pixel plans sizes sib 

Jd "V+Jh-I - 

K D = K S +K B -1 

Hnrimftd impulse response array node: 

The 2D impute response function array is rotated 180'tegreo and translated with.respect to the j-k 
pixel plane of the source image. A destination pixel is computed whenever all elements cf Hie array lie 
within (he boundVof the source image; all other •dwttoaticm'pbtels'. are set to rero.'-'Tfce* destination 
image size is the same as the source image size. 

For all 1. xn, n. and for ' 

j« i j £ if J H .isodd . 

Jc-* £J £. J-l-j e if J H is even 

i k £ K^l^k, ifKaisodd 

k, -1 < k < X-l-k, if K R is even 



DO".k.I.m,n) , 



P £S (j', 1. a.n)H{j - j* +j c . k - k' + k c ) 



For all 1. m, n, end for . ... . " "• " 

j < J e ' -••ar'j- > J-l-j c - ^ ■- - •if^isodd. . 

j<j 0 -l or j >J-l_j e ifJ H iseven 

k< k. or k > K-l-k, ifKaisodd 

k<k e -l ork>K r l-k c ifKHiseven 

DQ, k. 1. in. nH 0 . 
where ■ ■ • 



j 



J H -1 



. ,. . ifJ«lsodd 

J — n 
t ~ ~2~ if Jjj is even 
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•-• convolveJZd 

*c a —5 — m . if K« is odd 

*c *" ~~T if K H is even 

The summation limits fire 

MAJC{0.j-j o ) * j" S MlN{W.j+j a ) ifJ„isodd 

MAX{0.j-j B+ l}i j- i MEN{J-l,j+j e } .. if J B ii even .. . 

. MAX{0. k-JCc) £ k*S lUNIK-l.***^ if K s is odd ,• . 

" . \ ..^AXib.^+lj s _ k^ -.MI^iK-l.k+k,} ,' J ifKHise^en; ' 
He pfceel plate sizes are 

Kj> =*Kg — • K 
= *L.$ = L 

N D - N s ^ N 

Key pixel, zero exterior mode: 

^^,r!?£ "T"! fuaCdC ? f^ y iS r£,tated J8 ° dc ^ «* Elated with aspect to the j-k 
pLtcl plane of tic source image A destination pixel « compute whenever the k*y. pixel rf the array b 
over a source pixel. The destawtioa imaae n** U ti„. e «,JT» .i.- y K 1 _ arra J' » 

— • -j .t_ ■ , , . wnage size is lae same as ue SOurca image size Source mxrfc 

outside the j-k plane boundary ere assumed iq be zero valued P 

For all 1, m, a. and for 

MAX[0. < j ^ MIN{J-l,J~l+y 
MAXiO.-IqJs k i MtNlK-l.Jk-l+kj.) 

D(j ( k.l,m,n) = £ pSG 1 .kM.m,n)H(j-'j-+j v k-.k' + i: k ) 
For all 1, m. a. oad for 

j < MAX[6, -j,.'} or j > MtN{J-l.J-l+j t J 

k< MAXia-kt) or k > MU\ T (X-i;K~l+k t ) 

D(j. k. U. n) » 0 
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.convoiveJW 

The summation, limits are 

MAX(0.j-Jirtj fc +l} & y s MIN{J~l.j+j k } : , 
MAXtO.k-KH+kfe+l) ifc'i M3N{K-I,fc+J^}' 
The pixel place sizes are 

KjJ w Kg = K 

L D = Lj = L • . .. , 

N 0 - N s * N 



— — ^ - . .. /,'„ es 

N..|Cey pixel, reflected exicrigr modeL.-'' 
. — ' 

TTiB"2DTfflpuise response function array » routed 180 degrees and translated with Aspect -to the j-k 
pixel plane of the source image. A desdnadon pixel Is computed whenever the key pixel of the impulse 
response array is over a source pixel. u» destination image size is the same as the source imaresize 
Source pixels outside the j-k pixel plane boundary are assumed to be reflections of souxee pixels about 
each boundary, i.e., the source image is assumed to be reflected about its left, right, top. and bottom 

For all J, m. n. and for 

o < j < J- 1 
0 £ k £ K- 1 



33 (j. k.l,a. a) 

where 



j" "' 1 "" J* J' < 0 k" B 1 _ fc' k' < 0 

J" - J" k"- K OsV sx-i 

j- - (2CJ - « - j* j* > J - J k" * 2CK-])-k" k'>K-l 



Tha suTmnarion limits are 

The pixel plane ***** on* 

Jd - J S « J - 

Ld « I* = L 
Mi> 14 Ms = M 
Nx> ■ N s =1 N 

Profile: . „•'• 

, Foundation.. \; ; ^ " " ' '* '** * 

' - t • •. s ; . , ■ 

Error codes: 

•» # , 

001 
003 
Oil 
012 
013 
014 

041 • 
045 
052 
061 
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m& r/]\}fe s-j 
UiJKft rcQHJIKttJ 



*. i j -l^J aA, P*W* <*K hues X3- z.<* 

( tnfthsh o^t+^Ud hetetv) 



j 



MM: 



IP* 

life' 



1 
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Xo 

3(x!-xo) 
3(x a '-3xi+2x«) 



^•(2x 3 +7% +3x0 
^-(2x, + x a -3x,) 
-|"(2x 3 -5x a +3x l ) 



-^(xH^Xs+x*) 

-i-(X4-X I ) 

x 4 ~*2x3+Xj 



-|-(x 9 +4x<+Xj) 
a:,'— 2X<+Xj 



KI: *<(«) = V° + V <> « + V ) « a +6i co « 3 <i«1.2,-».N-1.0S«Sl) (15) 
till a, (») littM c^tllfc^i i^fiotttt. y 1***1 ti^^*c t> W#fc» * 

MB*. *r. 



x,(0) = W> = 
x',(0) « V° 



= -yC^f+i — X<-i) 



(16. a) 
(16.b) 

(16. c) 
(16.d) 



10: 

tit 

Mu) - j« J («-l)i w +|«(l+4«-3« l )J w +}(«-^ 3l<, - 2!( " 2)j: ' 



— ^* a (w-l)x,-i (i = 1.2,-.JV-l. 0S«S1) 



(17) 
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, s 'Hi 

*„(*) = V 0, +6i <0> » + W 0> " z+ *3 (0,w3 (-^"^ 0) (18) 

^(-D = v o> -V 0) +^ (fl> -*s (0) » * (19a) •'Ig 

* 0 (0) = b Q ™ - x, (19- W 

»'(0) - V 0) = T U_ ^ } (19 ' C) *!l 

= 6i<°>+26 a <°>+3V 0) - -yU-*i> (19. d). 

*(«> - ■^»(l + «)a*+^«<i*+l>CB-3k)ak+^(i*+l)(3ii 1 -10« + 10)« :J| 
1 u{u*-4u+5)x 0 <-l*«*0> (20) 'II 



A--.*. 



*"» - 'a 



xAu) = ^(«-l)(w^+2« + 2>^-^-(w-2)(3u 2 +4w + 3)^-i ^ 
+-i-(«-l) («-2) (3«+2)x*-z-~-U~l) i («-2)x*-3 (21) -p/ 

(1^M^2) 



,c.< 



List3.l.5^t*'f. |i]-/n^*9 set-blend 0 Wrt^€-lS^v»T. J 
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<List3.i.5> gMa^^7^'^^7/> (&^-^o<?<) 



4^1 . 



-ft!*. 



#include <stdio .h> 

# include <s tdl ib ♦ h> 

#include <conio.h> 

#include <graphics.h> . 

ftinclude "spline. h w 

#define SPLIT 8 
♦define X_RANGE 0.8 
#def ine Y_RANGE 0 ,8 

double xltlOOJ, yl[100], wx[SPLIT+l], wy[SPLIT+l), d; 
double b[SPLIT + l][4], ^^J"""** 1 ' ^ 

void set blend (double <*b>[4], double (*f>(4], double (*1)[41, int n) , 
void ink blend (double* x, double* y, double* wx, double* wy, 
double (*b)[4], int n); 

void main() 

double kx, ky, xmax, ymax, xmin, ymin; 
int xO, yO, y_height, x_width, n, k, cl; 
char *s~( "datal") ; 



*/ * 



n => rd data(s, xl, yl, sxmin, fixmax, 6 ymin, syntax)./ 

~ 7* read data from a file 1 

opengraph C'a: WtcWbgi." / *x_width, fiy_height) ; 

cl = getmaxcolor () ; 

kx x width * X_KANGE / (xroax - xmin) ; 
ky ~ y~height * Y_RANGE / (ymax - ymin); 
x0 =» x_width * ( 1-X_RANGE) / 2 - xmin * kx; 
y 0 o y_height * (l-Y_ftANGE) / 2 - ymin * ky/ 

plot_data(xl, yl f n, cl, xO, y0, kX/ ky, yjieight); 

draft spline sample program for open curve */ 

set blend (b, first b, last_b, SPLIT); /* set coeffcients */ 

k a T 0; { *~ /* data counter */ 

mk blend (&xl[k], syltk], wx, wy, first_b, SPLIT); 

~ /* first section */ 

mk,curve(wx, wy, SPLIT, cl, xO, yO, kx f ky, y_height) ; 

whilst k < n-3 ) t 

mk_blend(«xl[k] , &yl[k], wx, wy, b, SPLIT)*; 

rok_curve(wx, wy, SPLIT, cl, xO, yO, kx, ky, y_height) ; 

k++; 

} 

k— / 

mk_blend(&xl£k], &yl M # wx, wy, last Jd, spLXT) ; 
mk_curve(wx, wy, SPLIT, cl r xO, yO, kx # ky, y_heig*it); 
getch ( ) ; 
closegraph( ) ; 



void s etJ>lend (double (+b)[4], double (*f > [4] , double (*l)[4l, int n) 
{ 

double u; 
int i; 

fori i - 0; i <- iw i++ ) ( 
u » (double) i / (double) n; 
btiHO] - -u * (u - 1.0) * (u - 1.0) / 2.0; 
b I 1 - (U " 1.0) * (3.0- * u * u - 2-0 * u - 2.0) / 2.0; 
b[i][2] - u * (1.0 + 4.0 * U - 3.0 * u * u) / 2.0; 



bti] 13] - u * u 



(u - 1.0) / 2.0; 



u » (double) i / (double) n - 1.0; 
flillO] « -u * (u * u - 4.0 * u + 5.0) / 10.0; 
f (i] ID = (u + 1.0) * (3-0 * u * u 
f[il[2] - u * (u + 1.0) * (5.0 - 3.0 * u) / 10.0; 
f [i] [3) = u - u * <u + 1.0) / 10.0; 



- 10.0 * u + 10.0) / 10.0; 



u - (double) i / (double) n + 1.0; 

l[i](0) = -(U - 1.0) * (u - 1.0) * (u - 2.0) / 10.0; 
llilll] =» (u - 1.0) * (u - 2.0) * (3.0 * U + 2.0) / 10.0; 
1 i)I2J - -(« - 2.0) * (3.0 * u * u + 4.0 * u ♦ 3.0) / 10. 0; 
lti][31 - (u - 1.0) * (u * u + 2.0 * u + 2.0) / 10.0; . 




rD/7A*List3.I.BtwL*t. Wft set-blend 0 ifiMi^ List 3 .1 .5« *> 



<Llst 3 . 1 . 6> M L* £ ffltla) tz «> 0)S»£^ * 9 -f > • MtA 



/* draft spline sample program for closed curve */ 

set_blend{b, firstjb, lastjs, SPLIT); /* set coeffcienta */ 
^ a 0/ ~" /* data counter */ 

while( k < n-3 ) { 

mk blend ( sxl t k ) . syllkl, wx, wy, b, SPLIT) ; 

mk~curve(wx, wy, SPLIT, cl, xO, yO, kx, ky, yjvsight) ; 

k++; 

) 
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Chapter 3 APPLICATION PROGRAMMING 
3.1 Interpolation and Spline 



EXTRACTION 

<Table 3.1.2> Values of x, x' , and x" of starting four points 
of base spline 



n 


0 


. 1 


2 


3 


X 








-^■(2x3 + 7^^337,) 


+ 4x 3 +.r 2 ) 




x' 


3(Xi- 


-Xq) 






y(x,-x*) 


-g-(x 3 -x 3 ) 


x« 




-3x l +2a<,) 




-^(2x 3 — 5x a +3xi) 




X$ '— 2x* "1" x 3 



Curve rule-lxke spline 

When one dot and another are interconnected by using 
a curved rule,, they are interconnected in such a manner that 
an interconnecting line may be smooth at each of these points. 
To be "smooth" means that a linear differential coefficient at 
a data point is the same in the immediately preceding interval 
and in the immediately following interval. Of course, the line 
passes through the data point. Such a method of interpolation 
is referred to as "curved rule-like interpolation' 1 '" and as 
"curved rule-like spline" here. 
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Assume four consecutive data pieces to be X^, X x , X ± + 19 
and X t + 2 (where y-axial data is omitted) • Assuming a normalized 
parameter to be u(0£u£l), the following cubic equation is 
considered: 

zt(u) - V ) +6i <<, « + V°« 8 + W (i=1.2."'-,iV-1.0£K^l) (15) 

where X ± (u) is an x-axial value obtained by interpolation. A 
y-axial value can also be obtained similarly. First, to obtain 
coefficients b 0 U) , b^ , b 2 U) , and b 3 l » , conditions are given as 
follows : 



*,(0) = K« = Xt (16a) 

= V°+6. c « + V 0 + V = xt +l ' (16. b) 

= = -£<*«-*,..,> (I6.c) 

*',(!) = W+ZW+W* - (16. d) 



where linear differential coefficients at the respect points 
of Xi and X itl provide coefficients in a case where a quadratic 
equation is applied to three equally spaced points. By 
obtaining each of the coefficients from Equation (16) and 
organizing it, the following equation is obtained: 

(i-l t 2, : - t N-l. 0£tt£l) (17) 
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Equation (17) cannot be used for a first interval and 
a last interval* This is because u=-l and u=2 do not give points 
Xi_i and X im2 respectively (although these points have been passed 
though because a mixed function gives the same shape as that 
of Lagrange's interpolation). Therefore, for both ends, 
different conditions are given. That is, assuming an 
interpolation that gives a first interval to be X 0 (u) and 
performed in a range of -l£u£0, the following equation is 
given: 

Xa(u) = V w +V 0) ttH-V w » a +V 0, K a C-ISmSO) (18) 

The conditions are such that a starting point X 0 and a 
point X t should be passed through and a linear differential 
coefficient at the point x t should agree with that of an 
adjacent interval. Since these conditions alone are not enough 
to determine a coefficient, another condition is set such that 
the differential coefficients should agree also at a point X 2 . 
The conditions are written as follows: 



^(-1), = V 0) -&i <0> + £ 2 (0) -fe< 0 > 
*o(0) = V 0) 

xb'(l) = 6, (<>) + 2& ^ <0, ^-3V 0, 
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(19. a) 
(19. b) 

(19. c) 
(19. d) 



These are solved and organized in terms of x as follows: 
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Xo («) = -^-k 2 (1 + u) x 3 +±u (u+1) (5- 3k) & +Jq(u + 1) (3 w 2 - 1Qk+ 10) 
— ^w(« i, -4i<+5)ab (-l^w£0) (20) 



Similarly, in an ending interval, the following is 
obtained: 



x*(k)- ~ (w 3 +2w-l-2)a:w-^-<«-2) (3« 2 +4« + 3)x^_ 



+jq(u-1) (w-2) (3w + 2)x*-2~~)(w-l) a (i<-- 2)^-3 (21) 

(1^«^2) 

This interpolation by use of a spline gives results 
similar to those by a mixed function. A relevant program also 
takes on a similar form. Although the mixed function may seem 
to be better just because the equation is simplex in form, it 
may be common knowledge that the mixed function has such 
properties that the quadratic differential coefficients agree 
and it lacks intuitiveness, whereas by this spline the linear 
differential coefficients agree at data points. 

< 

• Programming of curved rule-like spline 

Programming of curved rule-like spline is indicated in List 
3.1.5. This programming is the same as that of a mixed spline 
except contents of a function setjblend() • Special processing 
is required only in a starting interval and an ending interval. 
The function mk_blend() is the same and so omitted in the list. 
Based on data which has been used as an example, almost the same 
curve as that by the mixed spline is obtained. 
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<List 3,1 ,5> Curved rule-like spline program 

(continued to next pag© ) 



# include <stdio.h> 

# include <stdlib.h> 

# include <conio.h> 

#include <graphics.h> . 

# Include * spline . h " 

# define SPLIT 8 
#define- X_RANGE 0.8 
#define Y~RANGE 0;8 

double xl[100), yl[100]', wx[SPLIT+!J, wy[$PLXT+l] , d; 
double btSPLlT+1] (4] , f irst>(SPLiT+l] [41 , . last_b[SPLtT+l] [4] ; 
void set_blend (double <*b)[4], double (*f)[4], double <*1)[4], int n) / 
void mk_blend (double* x, double* y, double* wx, double* wy, 
double (*b)[4], int n) ; 

void oiainl) 

double kx, ky, xmax, ymax, xrain, ymin; 

int xO, yO, y_height, x_width, n, k, cl; 

char *s«{"datal") ; 

n » rd_data(s, xl, yl, sxmin, &xmax, &ymin, «ymax); 

/* read data from a file 's' */ * 
opengraphC'aiWtcWbgi", &x_width, $y_height); 

cl a getmaxcolor ( ) ; 

kx » x_width * X_RANGS / (xaiax - xntin) ; 

ky « y_height * Y_RANGE / (ymax - ymin); 

xO = xjwidth * (1-X_RANGE) / 2 - xmin * kx; 

yO « y_height * (1-Y_RANGE) / 2 - ymin * ky; 

plot_data(xl, yl, n, cl, xO, yO, kx, ky, y_height) ; 

/* draft spline sample program for open curve */ 

set_blend(b, first_b, lastja, SPLIT); /* set coeffcients */ 
* = 0; 1 " /* data counter */ 

mkJfc>lcnd<sxUk] , &yl[k], wx, wy, first_b, SPLIT); 

" /* first section */ 
mk_curve(wx, wy, .SPLIT, cl, xO, yO, kx, ky, y_height) ; 
while ( k < n-3 ) { 

mk_blend(*xl [k] , «yl[k], wx, wy, b, SPLIT); 

mk_curve(wx, wy, SPLIT, cl, xO, yO, kx, ky, y_height) ; 

k++; 

) 

k~; 

mkjDlend(*xl [kj , syl[k), wx, wy, last_b, SPLIT); 
rak_curve(wx, wy, SPLIT, cl, xO, yO, kx, ky, y_height) ; 
getch(); 
closegraph() * 
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void setjplend (double <*b)[4], double (*f)[4], double (*1)[4], int n) 
< 

double u; 
int ±7 

for( i « 0; i <= n; L++ ) { 
u = (double )i / (double) n; 

b[i] [0] = -u * (u ^ 1.0) * (u - 1.0) / 2.0; 

[1] = (u - 1.0) * (3.0 * u * u - 2.0 * u - 2.0) / 2.0; 
b[i] [2] - u * (1.0 + 4.0 * u - 3.0 * u * u) / 2.0; 
t>[±] [3] « u * u * (u - 1.0) / 2.0; 

u = (double)! / (double )n - 1.0; 
f (±J [0] o -u * (u * u - 4.0 * u + 5.0) / 10.0; 

ftillD = (u + 1.0) * (3.0 * -u * u - 10.0 * u + 10.0) / 10.0; 
f(±J[2] - u *■ (u + 1.0) * (5.0 - 3.0 * u) / 10.0; 
£14] [3] = u » u * (u + 1.0). / 10.0; 

u « (double )i / (double >n +1.0; 

1[4][0] - -*(u - 1.0) * (u - 1.0) * (u - 2.0) / 10.0; 
l(i][l] - (u - 1.0) * (u - 2.0) * (3.0 * u + 2.0) / 10.0; 
1UJC2] - ~(u - 2.0) * (3.0 *u*u + 4.0*u + 3.0) / 10. 0; 
114} [3] ~ (u - 1.0) * (u * u + 2.0 * u + 2.0) / 10. 0; 



} 

) 
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• Program for closed curve 

A program is indicated in List 3.1.6- It is the same 
as that of the mixed function except that the function 
set_blend() is such as given in List 3.1.5. 

<L±st 3.1. 6> Curved rule-lijce spline program for closed curve 

/* draft spline sample program for closed curve */ ■ 

set_blend<b, first_b, last_b, SPLIT); . /* set coeffcients */ 



3c - 0; 

while ( k < ti-3 ) { 



/* data counter */ 



»kJ&lend<fixX[k] , *yl[kj, wx, wy, b, 
mk_cuxve(wx, wy, SPLIT, el, xO, yO, 
3c++; 



SPLIT) ; 

ta ' ty' y_height); 



getchO ; 
closograph () ; 





• Execution example* 
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